Orientations of two coupled molecules 
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Orientation states of two coupled polar molecules controlled by laser pulses are studied theoreti- 
cally. By varying the period of a series of periodically applied laser pulse, transition from regular to 
chaotic behavior may occur. Schmidt decomposition is used to measure the degree of entanglement. 
It is found that the entanglement can be enhanced by increasing the strength of laser pulse. 
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The ability to control the alignment of molecules may 
have important purpose in stereodynamics surface 
catalysis trapping molecules 0, molecular focusing 
01 , and nanoscale design |^ @ . Experimentally, several 
methods have been used to control the orientations of 
molecules 0,11, 0,0, ^2- For example, by turning on a 
picosecond laser pulse adiabatically, the pendular states- 
hybrid of field-free molecular eigenstates 0, 0, Q| - 
can be created. A femtosecond laser pulse, like impulsive 
excitation, is found to be able to generate a field-free 
orientation T?,"]^. 

Since entangled states are fundamental for quantum in- 
formation processing 0, 0| , many research works have 
been proposed to generate entanglement in quantum- 
optic and atomic systems . In this Letter we pro- 
pose a novel way to generate entanglement between two 
coupled identical polar molecules separated in a distance 
of tens of nanometers. Both molecules are irradiated by 
ultra-short pulses of laser light. Our study shows the 
entanglement induced by the dipole interaction can be 
enhanced by controlling the inter-molecule distance and 
the field strength of laser pulse. 

Consider now two identical polar molecules (separated 
by a distance of i? ). There exists dipole-dipole inter- 
action between these two molecules. Ultrashort half- 
cycle laser pulses are applied to both molecules . The 
Hamiltonian of the system can be written as 

where L| and |j (= B) arc the angular momentum 
operator and rotational constant, respectively. Udip = 
[Mi ■ M2 ~ 3 {jli • gfl) (^2 ■ ^r)]IR^ is the dipole interac- 
tion between two molecules, where Jli and JI2 are the 
dipole moments. For simplicity, the dipole moments of 
two molecules are assumed to be identical, i.e. = 
^2 = M- The field-molecule coupling can be expressed as 
Hi = —fiE {t) cos 9 cos (ojt) — (t) cos 6 cos (ujt) , where 
6 and 9' are the angles between the dipole moments and 
laser field. The laser profile is assumed to be Gaussian 
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shape, i.e. E (t) = Eqc ^ , where Eq is the field 
strength, to is the center of peak, and a is the pulse 
duration. With these assumptions, the time-dependent 
Schrodinger equation can be solved by expanding the 
wave function in terms of a series of field-free spherical 



harmonic functions 

^ ci^i'm'{t)Yi„,{9,(f>)Yi,,n,{9',(f>'), (2) 
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where {9,<j)) and [9',(j)') are the coordinates of first and 
second molecule respectively, cimi'm' {t) are the time- 
dependent coefficients corresponding to the quantum 
numbers {l,m;l' ,m') and can be determined by solving 
the Schrodinger equations numerically. In above equa- 
tion, total wavefunction has no spatial dependence since 
we keep inter-molecule separation i? as a fixed parame- 
ter. One might argue that variation of R is inevitable 
because of the effects of laser fields or inter-molecule vi- 
brations. However, recent experiments have shown that 
it is possible to resolve two individual molecules sep- 
arated by tens of nanometers when they are hindered 
on a surface p^. |2^. In principle, the free orientation 
model can be easily generalized to hindered ones by re- 
placing the spherical harmonic functions with hindered 
wavefunctions as shown in our previous work 24]. The 
essential physics discussed below should be similar. 

The orientations (cos 6*) and (cos 6*') can be evaluated 
immediately after the coefficients cimi'm' {t) are deter- 
mined. The parameters for numerical calculations are 
based on Nal molecule whose dipole moment is 9.2 de- 
byes and rotational constant is 0.12 cm~^ in t he g round 
state. Other parameters are referred to Ref. hj. The 
duration and frequency of the half-cycle pulse are set to 
279 fs and 30 cm~^, respectively. The center of peak 
is 1200 fs and the initial condition is cqooo = 0) = 1- 
Fig. 1 shows the orientations of the first and second 
molecules after a single laser pulse is applied on both 
molecules. For R — Zx 10^^ m, the behavior of the first 
molecule is quite close to that of a free rotor [l^ . This is 
not surprising because the dipole interaction is weak for 
this molecule separation. However, as two molecules get 
close enough (Fig. 1 (b)), both molecules orient disor- 
derly, and the periodic behavior disappears. This is be- 
cause the dipole interaction is increased as the distance 
between the molecules is decreased, and the energy ex- 
change between two molecules becomes more frequently. 
The regular orientation caused by the laser pulse is in- 
hibited by the mutual interaction. 

The populations of some low-energy levels are shown 
in the lower panels of Fig. 1(a) and (b). The solid, 
dashed, and dotted lines represent the populations of 
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the states (1,0;0,0), (1,0;1,0), and (2,0;1,0), respectively. 
These states show different degrees of periodic behavior 
at different distances. However, the populations of some 
higher excited states, for example the (3,0;1,0) state in 
the inset of Fig. 1(b), display different degrees of irreg- 
ularity. This manifests a fact that the nonlinear effect, 
caused by the reduction of i?, does not affect the regu- 
larity of the low- lying states, and the origin of the irreg- 
ularity is caused by the higher excited states. 
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FIG. 1: Upper panels of Fig. 1(a) and (b) show the ori- 
entations of the two molecules at different distances. Lower 
panels: The populations of the states {I, m; I' , m') — {l, 0; 0, 0) 
(solid lines), (2,0; 1,0) (dotted lines), (1,0; 1,0) (dashed 
lines). The insets in (a) and (b) represent the population 
of state (3,0; 1,0). 

Consider now the molecules are irradiated by a series 
of laser pulses periodically. As shown in Fig. 2 (a), if the 
period of the applied periodically laser pulse T is equal 
to Ti/B , then both molecules behave disorderly no mat- 
ter how the distance R is varied. The chaotic behavior of 
the molecules can be ascribed to the well-known "kicked- 
rotor" problem. However, a series of regular-like orienta- 
tions marked by dotted and dashed lines are present in 
Fig. 2 (b) if T is equal to ttTi/B. For a free rotor under a 
single kick, this interesting phenomenon comes from the 
situation as the magnitude of the orientation returns to 
its initial condition {{cos 9) — 0) after a certain period 
T [T^ . Therefore, for two molecules in weak interaction 
limit {R = 3 X 10~^ m), the wavepacket-like orienta- 



tion is similar to that of a single free rotor under the 
same laser period. The difference is the suppression of 
the amplitudes at long time (dashed lines). It means that 
the dipole force can generate some accidental phases to 
perturb the regularity of the coupled system. The lower 
panel of Fig 2 (b) exhibits that the suppression of the 
regularity is quicker if the dipole force is stronger. 
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FIG. 2: The orientations of the first and second molecules 
under periodic laser pulses with the intensity Eg — 3 x 10^ 
V/m and periods T= (a) Ih/B, (b) nh/B ps. The upper and 
lower panels of (a) and (b) correspond to the distances R — 3 
X 10"* and 2 x 10~* m, respectively. 

Let us now turn our attention to the entanglement gen- 
erated in this system. The coupled molecules in fact can 
be expressed as a pure bipartite system 



(3) 



Iral'' 



The partial density operator for the first molecule is given 

by 



Pmoi 1 = Tr„,oi 2 I*) (^-1 



(4) 



Following the procedure of Schmidt decomposition, the 
bases of molecule 1 is rotated to make the reduced density 
matrix Pjjjoi i diagonal. The entangled state can then be 
represented by a biorthogonal expression with positive 
real coefficients as given by 



I*) 



l-^™)mol 1 \Ylrn) 



™/mol 2 ' 



(5) 



where A/„j is the eigenvalue corresponding to 
l^'™)raoi 1 l^'™)raoi 2- The mcasurc of entanglement 
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for the coupled molecules can be parametrized by the 
von Neumann entropy 



Entropy = - ^ Ai„log„A;„ 



(6) 



Im 



Fig. 3 shows the time-dependent entropy after one pulse 
passes through this system. For inter-distance R = 5 x 
10~^ m, the entropy increases slowly from zero. For B = 
1.5 X 10^^ m, on the contrary, the entropy grows rapidly 
with the increasing of time because the dipole force is 
stronger. Notes that the entropy only varies within a 
finite range at long time regime. This indicates that the 
system reaches a dynamic equilibrium state even though 
the dipole force is still present. 
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FIG. 3: Time evolution of the entropy after applying single 

laser pulse for (a) i? = 5 x 10"* m and (b) 7? = 1.5 x 10~* 
m. The field strength is set equal to 3 x 10^ V/m. 

Fig. 4 (a) illustrates the variations of the entropy with 
respect to different field strengths of the applied laser 
pulse as R is set equal to 1.5 x 10~^ m. For the field 
strength Eq = 1.5 x lO'^ V/m, an irregular- like behav- 
ior of the entropy is obtained, and its value is not large 
enough for quantum information processing. However, 
Fig. 4 (b) shows that the degree of entanglement can 
be enhanced if one increases the field strength. This can 
be understood well by studying the relationship between 
the dipole interaction and the field strength. If the effect 
of dipole interaction overwhelms the laser field, most of 
the populations are distributed on the low-lying states. 
In this case, the entropy from Schmidt decomposition is 



certainly small as shown in Fig. 4 (a). On the other 
hand, if the field strength plays a dominant role, the dis- 
tribution of molecular states covers a wider range and 
the entropy is enhanced in this limit. 
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FIG. 4: Time evolution of the entropy for inter-molecule sep- 
aration R = 1.5 x 10~* m. The degree of entanglement can 
be enhanced if one increases the field strength. 

A few remarks about the differences between present 
proposal and previous works on generating entangled 
states should be emphasized here. In our model, the rota- 
tional excited states instead of internal electronic states 
of the molecules are considered. Second, the laser fre- 
quency in our work is tuned far-away from resonance, 
while conventional creation of entanglement depends on 
the resonant driving pulses. This means our work pro- 
vides a wider range to select the laser frequency to create 
entanglement. As for the effect of decoherence, our en- 
tanglement is formed by the excited rotational states, 
instead of the vibrational states. Therefore, the deco- 
herence is dominated by photon emission even if the 
molecules are attached to the surface of a solid. 

In conclusion, we have studied a system of two coupled 
polar molecules, irradiated by laser pulse. The chaotic 
behavior of the orientations comes from the populations 
in the excited states of higher mode numbers. The entan- 
gled states vary with the field strength of laser pulse. And 
more stable entangled state can be generated even though 
the coupled system has irregular orientations. This fea- 
ture may be useful in quantum information processing. 
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